arXiv:hep-lat/9808035vl 25 Aug 1998 


i 


The QCD Dirac Operator Spectrum and Finite-Volume Scaling 

P.H. Damgaard a * 

a The Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 

Random matrix theory and chiral Lagrangians offer a convenient tool for the exact calculation of microscopic 
spectral correlators of the Dirac operator in a well-defined finite-volume scaling regime. 


1. INTRODUCTION 

The succesful analytical computations of the 
so-called microscopic QCD Dirac operator spec¬ 
trum from random matrix theory 01 have 
raised an important question. Why can random 
matrix theory, which seems to be very far from 
the quantum field theory of QCD, apparently be 
used to compute spectral properties of the Dirac 
operator? Within the last year this issue has been 
much clarified, with important consequences also 
for lattice gauge theory. The results obtained 
in random matrix theory are now seen to be ex¬ 
act predictions of QCD itself, in an appropriate 
finite-volume scaling regime. This is defined by 
I 1/A qcd « L « 1 /m„, where L = V 1//4 
is the linear extent of the four-volume V. In 
this region one is probing finite-volume effects of 
QCD to the extreme. Only static modes of the 
lowest hadronic excitations contribute to the eu¬ 
clidean partition function of QCD, and if chiral 
symmetry is spontaneously broken according to 
SU(Nf) x SU(Nf) -> SU(Nf), it reads @ 

Z v = f dU(detU) u e VSReTr[Muf] (1) 

in a sector of fixed topological charge v. Here 
£ = (’ipip) is the quark condensate, and M is the 
quark mass matrix. In this context it is most 
useful to view Z v as the generating function of 
i/jiIj (because this is how one sees that it contains 
information about the Dirac operator spectrum). 
It depends on the masses only in the particular 
combination /x = mV £. This makes Z v a finite- 
size scaling function. 
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1.1. The double-microscopic limit 

Eigenvalues X n of the Dirac operator are deter¬ 
mined by ilj(t' n — X n (\> n . By comparing (jlj) with 
the original path integral (after integrating out 
the fermions), 

=n m / / i dA \^ nv™+ m }) e ~ s[A] . ( 2 ) 

/ J f,n 

one infers immediately that also the non-zero 
eigenvalues A in this finite-volume limit can only 
enter in the combination £ = AUE. It is thus a 
double-microscopic limit A, to —► 0 as V —>oo with 
fj, and £ fixed that will correspond to this finite- 
volume scaling regime [|]|5]||. Remarkably, in this 
limit random matrix theory can be used to give 
am entirely equivalent description of the partition 
function ©■ The relevant random matrix theory 
is given by {§ 

Z„ = fdWY[det(M + m f ) e -f*r e(m 2 ) ( 3) 
J f 

where 

" - ( T ) <‘> 

and IT is a rectangular complex matrix of size 
Nx (7V+|i/|). This matrix ensemble, which is the 
one relevant for QCD, is known as the chUE. The 
choice of ensemble, and the determinantal struc¬ 
ture in (|j), is fixed by the symmetry properties 
of the Dirac operator Q, while nothing a priori 
fixes the matrix potential V(M 2 ). For simplic¬ 
ity, denote also the eigenvalues of M by A. In 
ref. [0 it was shown for the particular case of 
a Gaussian random matrix potential that in the 
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double-microscopic limit the two partition func¬ 
tions agree (up to an irrelevant mass-independent 
constant): 

Z v (fi lt ...,n Nf ) = Z v (m,...,n Nf ) . ( 5 ) 

The required correspondences between the two 
formulations are E <-> 7rp(0) and V <-> 27V, where 
p( A) is the density of the random matrix eigen¬ 
values. One particular consequence of a series of 
universality theorems (see ref. [Q and the first 
of ref. ||) is that this equality holds universally, 

1. e. independently of the choice of random ma¬ 
trix model potential ||. The only essential re¬ 
quirement is, as one could expect, that p{ 0)^0. 
This means that masses and eigenvalues can be 
rescaled to common values.^] 

The universal equivalence © is the root from 
which all exact relations between random ma¬ 
trix theory and QCD in the finite-volume scaling 
regime 1/A qcd « L « l/m^ follow. Because 
there are long-established techniques for comput¬ 
ing spectral correlation functions in random ma¬ 
trix theory, this equivalence is not just interesting 
from a conceptual point of view, but is of great 
practical value. 

2. DIRAC SPECTRA AND PARTITION 
FUNCTIONS 

One particularly simple way of computing the 
microscopic spectral correlators of the Dirac oper¬ 
ator eigenvalues is by relations to partition func¬ 
tions with additional quark species. In random 
matrix theory based on the chUE ensemble there 
exists a convenient shortcut for the calculation 
of spectral correlation functions which is based 
on the kernel K N {z, z'\ mi,..., TOjv,). From this 
kernel alone one can derive all spectral correlation 
functions in the limit N —> oo: 

p(Ai,..., \ n , m \,, mN^ ) 

= det K(X a , \b]mi,... ,m Nf ) . (6) 

a,b 

This relation is exact, and not restricted to the 
double-microscopic limit. When comparing the 

2 It is tempting to consider the case of p(0) —> 0 in random 
matrix language as the limit of chiral symmetry restau- 
ration, but there is presently no obvious classification of 
quantum field theories that match the new random matrix 
theory universality classes with p(0) = 0 [fj. 


explicit expression for the kernel with the eigen¬ 
value description of the random matrix partition 
function ([|) itself, one finds a very simple re¬ 
lation involving the ratio of two partition func¬ 
tions, one corresponding to having two additional 
quark species of imaginary mass! In the double- 
microscopic limit this relation simplifies even fur¬ 
ther, and if we now make use of the equivalence 
©, we can express this kernel entirely in terms 
of finite-volume partition functions in the scaling 
region 1 /Aqcd << L «l/m n fb|: 

Ks((, C'; Mi, • • •, HN f ) =n yj (C 2 +M/)(C' 2 +M/) 
x^C, (7) 

A/(.Ml, . . . , PNf) 

up to an overall normalization constant that can 
be fixed by matching with the spectral density at 
the origin. Eq. (Q) is a master formula from which 
all microscopic spectral properties of the QCD 
Dirac operator can be derived. In particular, the 
microscopic spectral density itself is obtained by 
taking the coincident limit C = C /: 


ps{C ; pi,...,ii Nf ) = \C\ n(C 2 + M/) 

/ 

Z v (p 1 ,...,p, N iC,iQ) 

x —FT- L -f—’ ( 8 ) 

again up to an overall constant, which is easily 
fixed [jio]] . Although it involves a theory with two 
more quarks, one can also view it as the expecta¬ 
tion value of a complicated operator in the origi¬ 
nal theory with just Nf quarks. 

It turns out that also the distribution 
P.s((> Ai, • • ■, Av,) of just the smallest Dirac op¬ 
erator eigenvalue can be expressed in terms of the 
finite-volume partition function. In the sector of 
topological charge v — 0 it reads |]TT| 


PsiC, 

d_ 

"•dC 


Mi, 


-C 2 /4 


, P>Nf ) 


ZoW Mi+C 2 , • • •, yJl^Nf+C 2 ) 


Zc\( LL i 


11 AT 1 


( 9 ) 


and similar expressions can be derived for any v. 

The expressions (||) and ([jj) for the microscopic 
spectral density and smallest eigenvalue distribu- 
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tions are the most convenient quantities to mea¬ 
sure in lattice gauge theory. Their analytical ex¬ 
pressions, using the explicit form of the finite- 
volume QCD partition function [Q, agree pre¬ 
cisely with results obtained earlier based directly 
on random matrix theory computations 
Since the analytical results are known for arbi¬ 
trary masses (keeping the condition 1 /A qcd« 
L « they can in principle be measured 

by Monte Carlo techniques without difficulty even 
with dynamical fermions. Existing Monte Carlo 
meaurements ]l3| for QCD in 4 dimensions are 
based on staggered fermions and an SU ( 2 ) gauge 
group, for which predictions actually are slightly 
different. Very recently the expression (||) has 
been derived directly from the chiral Lagrangian 
framework, without relying on the relation to ran¬ 
dom matrix theory at all |l4|] . 

2.1. QCD in 3 dimensions 

One can also entertain the idea of spontaneous 
breaking of “chiral” symmetries in 3-dimensional 
QCD. Because fermions are described by two- 
spinors, and 7 -matrices thus can be taken as the 
Pauli matrices (and minus these matrices, a dif¬ 
ferent representation of the Clifford algebra here), 
there is no room for an analogue of 75 . How¬ 
ever, by pairing up quark fields one can refor¬ 
mulate the theory in terms of four-spinors. Two 
hidden chiral symmetries (based on 74 and 75 ro¬ 
tations) then become manifest. For an even num¬ 
ber of flavors it is conceivable that flavor symme¬ 
try, in the guise of chiral symmetry, can break 
spontaneously according to U(Nf) —*U(Nf/ 2) x 
U(Nf/ 2), a case that has been analyzed in the 
random matrix framework in p5[ . It recently 
been studied on the lattice in a quenched sim¬ 
ulation, and indeed a condensate was observed 
at zero temperature, disappearing at the decon¬ 
finement phase transition [ |l 6 | . In the broken 
phase one can study the microscopic Dirac oper¬ 
ator spectrum, which turns out to be somewhat 
different from the 4-dimensional counterpart. 

The case of quenched 3-dimensional QCD is 
the only known case for which the prediction for 
the microscopic spectral density is exactly flat. 
New Monte Carlo results on lattices of sizes rang¬ 
ing from 6 3 to 14 3 , show fairly good agreement 


with this prediction |16|| . Getting a flat micro¬ 
scopic distribution is in fact quite non-trivial, 
since generically the spectral density reflects the 
underlying peaks corresponding to the averaging 
over individual peaks of the smallest eigenvalues. 
It will of course be more interesting to see the 
result of an analogous simulation with dynamical 
fermions, for which the analytical prediction is a 
highly non-trivial oscillatory function |L5| . 
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